clc;
clf;
clear;
close all;

oam = lightSource();
oam.getTransOnceScale(1e-3, 512) % 只转换一次 + ，参数：半径，接收率
oam.getTransOnceScale(3e-3, 512)
oam.getTransOnceScale(5e-3, 512)
oam.getTransTwiceScale(1e-3, 512) % 转换两次次 +- ，参数：半径，接收率
oam.getTransTwiceScale(3e-3, 512)
oam.getTransTwiceScale(5e-3, 512)
% disp(oam);
oam.showLightIntensityDiagram(9e-3, 512);
% oam.showPhaseDiagram(9e-3, 512);
% atm = atmosphere();
% disp(atm);
% sim = simulator(oam, atm);
% disp(sim)
% sim.run();
% % sim.showLightIntensityDiagram(1000, 2, 512);
% sim.showLightIntensityDiagram(1000, 9e-3, 512);
% % sim.showPhaseDiagram(100, 9e-3, 512);
% sim.statisticalPhaseCrosstalk(1000, 9e-3);